Forced patterns near a Turing-Hopf bifurcation 
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We study time-periodic forcing of spatially-extended patterns near a Turing-Hopf bifurcation 
point. A symmetry-based normal form analysis yields several predictions, including that (i) weak 
forcing near the intrinsic Hopf frequency enhances or suppresses the Turing amplitude by an amount 
that scales quadratically with the forcing strength, and (ii) the strongest effect is seen for forcing 
that is detuned from the Hopf frequency. To apply our results to specific models, we perform a 
perturbation analysis on general two-component reaction-diffusion systems, which reveals whether 
the forcing suppresses or enhances the spatial pattern. For the suppressing case, our results explain 
features of previous experiments on the CDIMA chemical reaction. However, we also find examples 
of the enhancing case, which has not yet been observed in experiment. Numerical simulations verify 
the predicted dependence on the forcing parameters. 
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Turing patterns, originally conjectured as the basis for 



ally i 

biological morphogenesis [if, arise in such diverse fields 
as neuroscience 0] , ecology , materials science Q , and 
astrophysics Turing patterns were first observed 

experimentally in the CIMA chemical reaction [||; fur- 
ther experiments revealed a rich variety of patterns in- 
cluding spots, stripes, rhombs, hexagons, and supcrlat- 
tices 0, [1 ©, 03 . CIMA, and the closely-related CDIMA 
reaction, serve as convenient prototypes for other Turing 
systems. Additionally, these reactions are testbeds for 
questions of control. Spatial, temporal, and spatiotem- 
poral forcing have been shown to induce a transition be- 
tween patterns ll,|l2j, to introduce new (possibly local- 
ized) states 13, 14 1, or simply to suppress patterns fTHj ] . 



Many systems that form Turing patterns can also dis- 
play a Hopf bifurcation to time-periodic solutions. Fluid, 
chemical, and electrical experiments Til 17 . 18j and the- 



oretical studies [19|, [20|, [2l|, [22| have documented the 



rich dynamical behavior occurring when Turing and Hopf 
modes interact. If a system possesses Turing and Hopf 
instabilities and is susceptible to external forcing, one 
may be able to manipulate and control those dynamics. 

Thus, in this Letter, we study time-periodic forcing 
of spatially-extended systems near a Turing-Hopf bifur- 
cation point. We give special attention to the forcing's 
effect on the Turing pattern. The forcing drives the Hopf 
mode, which in turn enhances or suppresses the Turing 
pattern modes to which it is coupled. Our results are 
three-fold. First, we study weakly-forced Turing-Hopf bi- 
furcations from a symmetry-based (model independent) 
perspective. We predict the dependence of the solutions 
on forcing amplitude and frequency and tie our results to 
experiments on the forced CDIMA system. Second, be- 
ginning with a generic two-component reaction-diffusion 
system, we use perturbation theory to calculate coeffi- 



cients in the forced Turing-Hopf normal form. Of spe- 
cial interest is the coefficient of the term that couples 
the Hopf mode to the Turing mode, whose sign dictates 
whether the forcing enhances or suppresses Turing pat- 
terns. The result of the calculation can be applied to spe- 
cific two-component systems such as the Lengyel-Epstein 
model [23|, the Brusselator [24j|, and so forth to deter- 
mine, as a function of the reaction kinetics, the qualita- 
tive effect of the forcing. Though previous work on the 
CDIMA reaction and Leng yel-Epstein model found only 
the suppressing effect [la . |25| , we find enhancement in 
other models. Third, we perform numerical simulations 
which verify the symmetry and perturbation results. 

Symmetry analysis. — Consider a spatially-extended 
system with a spatially homogeneous base state (zero, 
without loss of generality) and a codimension-two point 
in parameter space where steady, spatially periodic (Tur- 
ing) modes and a spatially homogenous, time-periodic 
(Hopf) mode bifurcate simultaneously. For expository 
purposes, assume a simple stripe pattern. Close to the 
bifurcation point, let 



u = zr(^2)vTe lq ' : 



t H {t 2 )v H e WHt +C.C.+ 



(1) 



where u describes the system state (e.g., chemical con- 
centrations), ZT t ii(t , i) are the slowly-varying amplitudes 
of the Turing and Hopf modes v^h, t% is a slow time 
variable, q is the wave vector of the pattern with wave 
number q c = |q| as determined by linear stability the- 
ory, loh is the Hopf frequency, and the dots represent 
damped modes. Assume time-periodic forcing ./'(f). For 
sufficiently weak forcing, the only frequency component 
of f(t) entering the weakly nonlinear description will 
be that closest to u>h- Call this frequency ujf and its 
strength Define the detuning A = uif — loh- Wc 
assume A = for now; later, we relax this assumption. 



2 



Equations for the slow-time evolution of zt.h must re- 
spect the spatial symmetries of the underlying system, 
namely translation by Ax [zt — * zj'e !q ' Ax ) and inver- 
sion through the origin [zt — ► z~t)- Forcing breaks the 
temporal symmetry of the problem, but it may be re- 
cast as a parameter symmetry by incorporating the ac- 
tion of time translation by At on the forcing function 
(see, e.g., [1^1) so that the full symmetry is (zh, — ► 
( z h, I /h \)e iuJHAt . The cubic amplitude equations are 

Zt = Xz T + sign(g 1 )\z T \ 2 z T + g 2 \z H \ 2 z Tl (2a) 
zh = \fH\+^z H + gz\z H \ 2 ZH+g<i\zT\ 2 ZH- (2b) 

Here, A, gi^ € K and /i,g 3 ,4 G C. We perform a phase 
shift and rescaling so that the forcing appears as \fn\] 
we also scale such that the self-interaction coefficient of 
zt is sign(gi) = ±1. Throughout the symmetry analysis, 
we assume a supercritical bifurcation to stripes, and thus 
set sign(i?i) = — 1. For weak forcing, terms smaller than 
O(fn) have no effect at leading order and are ignored. 

In the absence of forcing (fn = 0) there are four so- 
lutions: the trivial solution, a pure Turing mode, a pure 
Hopf mode, and a mixed mode (see, e.g., [I(|). Weak 
forcing perturbs the solutions as follows. 

(i) Perturbed trivial solution. Here zt = 0, zh ~ 
— The physical solution |TJ is an O(fn) spatially 
homogenous oscillation with frequency ojh- 

(ii) Perturbed Hopf mode. In this case, 



zt = 0, zh ~ Re 



R(ti,g 3 ) = y/-Rc(fi)/Rc(g 3 ) 
n(n,g 3 ) = Im( M ) + Im( 53 )i? 2 



f H \(A + Be 2int ), (3a) 
(3b) 
(3c) 



The 0(1) constants A, B are obtained by substituting ([3]) 
into j2]). The physical solution is an 0(1) spatially ho- 
mogeneous oscillation with frequency ujh + O superposed 
on 0{Jh) oscillations with frequencies ojh and loh + 2fl 
The 0{Jh) correction is consistent with 271. l28l. |29| . 
(Hi) Perturbed Turing mode. We have 



zt 



ZH 



Vxh. 



92\f H \ 2 
2A|m| 2 

Jl = fi + \g 4 . 



(4a) 
(4b) 



The physical solution is an 0(1) spatial pattern super- 
posed on an O(fn) spatially homogenous oscillation. 
(iv) Perturbed mixed mode. In this case, zt and zh are 
both time-dependent: 



-H 



A + g 2 R 2 + \fn\C cos(ta 



Re 



'in 



\f H \ [A + Be 



2iS1t 



(5a) 
(5b) 



where R = R(fi,g 3 ),Q = fi(^,5 3 ) and g 3 = g 3 + g 2 gi. 
The 0(1) constants A, B, C, (j> are obtained by substi- 
tution. The physical solution is an 0(1) spatial pattern 



with 0{Jh) breathing at frequency fl, superposed on os- 
cillations similar to case (ii). 

To explore mechanisms for suppressing and enhancing 
Turing patterns, we consider further case (Hi). For com- 
parison with experimental results, we now allow detuning 
from ujh (A ^ 0), in which case (0| becomes 



zt 



Vx i 



ff2|/g| 2 

X\Ji-iA\ 



ZH 



\fs\e 



iAt 



iA 



(6) 



The forcing's effect hinges on the sign of g 2 - If g 2 < 
(g 2 > 0) the forcing suppresses (enhances) the Turing 
pattern. If g 2 < 0, weak forcing reduces the pattern am- 
plitude by a relative amount proportional to |/ir| 2 and 
shifts the domain of pattern existence from A > (un- 
forced case) to A > -# 2 |/iH 2 |M-£A|~ 2 > 0. If .92 > 0, the 
forcing enhances the pattern and shifts the bifurcation in 
the opposite direction, so that a pattern appears in the 
forced system for A values where it would not exist in 
the unforced system. The enhancement /suppression is 
strongest for A opt = Im(/I) ^ and decays away from 
this maximum. Note that since Jl depends linearly on fi 
and A, the shift of A opt away from is small if we are 
close to the codimension-two point. Furthermore, the 
rate of suppression (i.e., the eigenvalue of the perturbed 
Turing mode) is approximately — g 2 \fn | 2 A _1 \J1 — iA|~ 2 . 

The results above are consistent with experiments on 
the CDIMA chemical reaction [l5[ , which consider time- 
periodic forcing applied to stable Turing patterns. First, 
forcing was observed to suppress the pattern, correspond- 
ing to our case g 2 < 0. Second, suppression was strongest 
for ujf as ujh , in agreement with our explanation of cou- 
pling to the Hopf mode as the (indirect) control mecha- 
nism. Third, the rate of pattern suppression has a max- 
imum near A = 0. Finally, the numerical simulations 



of 15(, taken as a qualitative model of the experiment, 
showed that the domain of Turing pattern existence was 
shifted, with a critical curve in the (|//j|, A) plane having 
a minimum near A — 0, which follows directly from our 
expression for shifted domain existence. 

Perturbation analysis. — To connect the symmetry re- 
sults to models, consider a general two-component 
reaction-diffusion system with a spatially homogeneous 
steady state. As in [30j |. write the governing equations as 



R(u,«) + f(f), 



C 



dt~a-V 2 -b 

-c d t -d- KV 2 ' ' 



(7a) 
(7b) 



where C is the linear operator, f(t) is direct, spatially- 
homogeneous forcing and R contains nonlinear reac- 
tion terms. For convenience, define R 2 (R 3 ) to be the 
quadratic (cubic) terms in the Taylor expansion of R. 

From linear analysis, at the codimension-two point the 
parameters in (|7b[) satisfy d — —a, a 2 + be < 0, be < 0, 
a 2 (A' + l) 2 + AKbc = 0. The critical modes in Q arc 
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FIG. 1: Coefficients gi t 2 in <|2j as computed from two reaction- 
diffusion models, namely the Lengyel-Epstein equations (|9} 
and the Brusselator (|10|) . At the Turing-Hopf point, the co- 
efficients in ([2| depend only on the parameter A in the gov- 
erning equations. For each case, we focus on the interval of A 
for which g\ < so that the bifurcation to a Turing pattern 
is supercritical, (a) In the Lengyel-Epstein equations, gi < 
for A G [7.0, 13.8]. The coupling coefficient g2 < 0, so that 
forcing suppresses the spatial pattern, (b) In the Brusselator, 
gi < OforA € [0.68,3.26]. For A £ [0.68,2.57], g 2 < 0, so that 
forcing suppresses the spatial pattern. For A £ (2.57,3.26], 
(?2 > 0, so that forcing actually enhances the spatial pattern. 



v T = (-6, a - q 2 c ) T , v H = (~b, a 



T . The squared 
1)/(2K), and 



iujh 

critical wave number is |q c | 2 = q 2 = a(K 
the squared Hopf frequency is u) 2 H = —(a 2 + be) = Kq\. 

To compute the coefficients in ([2]), we perform a mul- 
tiple time scales perturbation expansion (for the case of 
weak forcing). The result of the perturbation expansion 
is (after rescaling) equation ([2]) , including expressions for 
the coefficients as functions of the parameters in (|7aj) . Of 
special interest is the sign of g 2 , 



sign p 2 = sign(v^, (<9R 2 £ 1 R 2 + R3 



(8) 



Here 



Ui is composed of the 

Y ■ S. 



4 = (c, -a + qif e^ 
critical modes in ((T|), and (r, s 

Consider first the Lengyel-E pste in model of the 
CDIMA reaction with forcing [Hl3l|, 



= A - u - Auv{l + u 2 )- 1 - f(t) 4 
= a [Blu-uvil + u 2 )' 1 + f(t)} 



VV 



(9a) 
DV 2 v\ (9b) 



where it, v represent the reacting chemical species and 
A, B, a, D are chemical parameters. At the codimension- 
two point, all coefficients in ^ can be written as func- 
tions of A. Fig. [TJa) shows gi,2(A). For all values of 
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FIG. 2: Simulation of the Brusselator (| 1 Of) . Values of it are 
indicated by shading in the x—t plane. The simulation begins 
with f(t) — (no forcing) and a random initial condition 
(not shown) which evolves to a steady state Turing pattern 
by t < 8000. At t = 8000, we set f(t) = 0.01 cos(Lu H t), which 
enhances the Turing pattern (as seen in the sharper contrast 
between peaks and valleys). The chemical parameters are 
A = 3, B = 9.998, D « 1.926 for which A = -Re/i = 0.001 
in ([2]). At the codimension-two point, q c « 1.47 and loh = 3. 



A for which the Turing pattern bifurcates supercritically 
(gi < 0), the coefficient g 2 < 0, and thus forcing sup- 
presses the pattern, in agreement with only suppression 
having been observed in CDIMA experiments. As a sec- 
ond example, consider the forced Brusselator f32| , 



u = A — (B + l)u 



= Bu — u v 



u 2 v 



DW 2 u, 



VV 



fit), 



(10a) 
(10b) 



for which all coefficients in ([2]) can again be written as 
functions of A. As shown in Fig. [l^b), when g\ < 0, 
there is one subinterval in which g 2 > 0, so that forcing 
will enhance the spatial Turing pattern. 

Numerical study. — Fig. [2] shows an example of Turing 
pattern enhancement in a simulation of (jlOjl . In further 
simulations we quantify the Turing pattern enhancement 
Et- To determine Et, we allow the forced system to 
settle onto its attractor. We measure the time-averaged 
amplitude of the critical Turing mode, and from this 
subtract the corresponding value for the unforced case. 
Fig. [3] shows Ex as a function of forcing strength for 
simple harmonic forcing and for three different sets of 
chemical parameters. The quadratic scaling of (O holds 
for sufficiently small forcing relative to the distance from 
the codimension-two point. Fig. [Ua) shows Et vs. A for 
the same three sets of chemical parameters; as expected, 
there exists an optimal detuning A op t that maximizes 
Et, and A opt — > as the system parameters approach 
the codimension-two point. Fig. BJb) verifies the func- 
tional dependence on A predicted by ([6]). 

In this Letter, we have studied the effect of time- 
periodic forcing on interacting Turing and Hopf insta- 
bilities. Our normal form results suggest a pattern con- 
trol mechanism for spatially-extended systems with these 
instabilities, including chemical reaction-diffusion sys- 
tems. Although the symmetry analysis is performed for 
weak forcing and small-amplitude patterns, it nonethe- 
less agrees with features of the experiments in 15jj . Fur- 
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FIG. 3: Scaling of Turing pattern enhancement Et as a func- 
tion of forcing strength for the Brusselator (|10|) with forcing 
f(t) = F cos(ojffi). Symbols correspond to simulations while 
the lines of slope 2 show the quadratic dependence on forcing 
strength predicted by (j6]). For all data, A = 3 and hence 
the critical wave number q c and Hopf frequency u>h are as in 
Fig. For each data set, A = -Re/i in ©. (•) B = 9.98, 
D w 1.935 for which A = 0.01. (■) B = 9.94, D w 1.955 for 
which A = 0.03. (A) B = 9.8, D » 2.031 for which A = 0.1. 




FIG. 4: Dependence of Turing pattern enhancement Et on 
frequency detuning A for the Brusselator (|10[) with forcing 
f(t) = F cos ((ujh + A)t). (a) The three sets of data corre- 
spond (same symbols) to the parameters in Fig. [3] except 
for the values of F. (•) F = 0.001. (■) F = 0.003. (a) 
F = 0.01. The lines guide the eye. As predicted by (|6]), there 
is a A opt for which is effect is maximal and A op t —> close 
to the codimension-two point, (b) Same data as • in (a), but 
here the line is a fit of the functional form ci/(c2 + (A — C3) 2 ) 
predicted by ©. 



thcrmore, we predict that the forcing may result in spa- 
tial pattern enhancement; thus far, only suppression 
has been observed in experiment. Nonetheless, we have 
demonstrated pattern enhancement with our perturba- 
tion analysis and numerical simulations. We hope that 
experimentalists might apply our results to look for the 
enhancing effect in other spatially-extended systems with 
Turing and Hopf instabilities. 
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